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An  Analysis  of  Combined  Longitudinal  and 
Torsional  Plastic  Waves  in  a  Thin-Walled  Tube1 

by 

Ro  J.  Clifton2 

Abstract .  This  paper  presents  a  one -dimensional  rate  independent 
theory  of  combined  longitudinal  and  torsional  plastic  wave  propa¬ 
gation  in  a  thin-walled  cylindrical  tube.  The  tube  material  is 
assumed  to  behave  as  an  isotropic  work-hardening,  elastic-plastic 
solid,  for  which  the  stress-strain  curve  in  simple  tension  is  a 
smooth  curve,  concave  toward  the  strain  axis.  The  resulting 
equations  are  shown  to  yield  two  wave  speeds,  c.  and  c  ,  which 
satisfy  the  inequalities  C£  -  cf  -  CQ  anc^  cs  -  c2  w^ere  c2 
elastic  shear  wave  velocity  and  cq  is  the  elastic  bar  velocity* 

The  solution  is  given  for  combined  longitudinal  and  torsional 
step  loading  of'  a  semi-infinite  tube,  which  is  initially  at  rest 
and  either  unstressed  or  statically  pre-stressed  to  arbitrary 
values  of  normal  stress  and  shear  stress.  The  solution  consists 
of  adjoining  centered  simple  wave  solutions  and  constant  state 
solutions.  There  are  two  types  of  simple  wave  solutions  corres¬ 
ponding  to  the  two  wave  speeds  c  and  c  .  The  stress  paths i in 

JT  S 

stress-space  associated  with  these  two  types  of  simple  waves 
form  an  orthogonal  network „ 

1  This  research  was  partly  supported  by  the  Department  of  Defense 
Advanced  Research  Projects  Agency,  Order  No.  71,  through  Contract 
DA-19-020-AMC-0077(R)  with  Ballistic  Research  Laboratories,' 

Aberdeen  Proving  Ground,  Maryland.  Partial  support  was  also  pro¬ 
vided  by  the  Advanced  Research  Projects  Agency  under  Contract  SD-86. 

2  Assistant  Professor  of  Engineering  (Research),  Brown  University, 
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Introduction 

A  test  involving  combined  dynamic  stress  states  which  is  attractive  from 
the  point  of  view  of  analysis  and  which  would  not  seem  to  present  unsurmountable 
experimental  difficulties  is  that  of  combined  longitudinal  and  torsional  impact 
of  a  thin-walled  cylindrical  tube.  Such  a  test,  for  stress  states  beyond  the 
elastic  limit  of  the  material,  could  be  expected  to  provide  useful  information 
regarding  the  dynamic  plastic  behavior  of  solids  under  combined  stress  states. 
Clearly  this  information  Is  of  fundamental  importance  in  the  development  of  a 
three-dimensional  theory  of  plastic  wave  propagation.  In  order  to  interpret  the 
results  of  plastic  wave  tests,  designed  to  study  the  dynamic  mechanical  behavior 
of  a  material,  it  is  necessary  to  first  solve  the  associated  boundary  value  pro¬ 
blem  for  an  assumed  material  behavior;  the  appropriateness  of  the  assumptions 
regarding  the  material  behavior  is  then  evaluated  by  comparison  of  theory  and  ex¬ 
periment.  Accordingly,  the  purpose  of  this  paper  is  to  present  an  analysis  of 
combined  longitudinal  and  torsional  impact  of  a  thin-walled  cylindrical  tube, 
based  on  an  assumed  material  behavior  which,  although  simple,  could  be  expected  to 
result  in  a  useful  theory  for  Interpreting  experimental  observations.  It  will  be 
assumed  that  the  stress-strain  cupve  for  simple  tension  is  smooth,  concave  toward 
the  strain  axis,  and  independent  of  the  rate  of  strain;  that  the  material  satisfies 
either  the  Tresca  or  von  Mises  yield  conditions;  that  isotropic  work-hardening 
applies;  that  the  elastic  and  plastic  strain  rates  are  separable  and  that  the 
plastic  strain  rate  Is  given  by  the  theory  of  the  plastic  potential.  In  other 
words,  the  theory  is  the  same  as  the  usual  quasi-static,  elastic-plastic  theory 
for  isotropic  work-hardening  (see  e.g.  Cl])  except  that  the  inertia  terms  are 
included o 

The  theory  presented  here  reduces  to  that  of  von  Karman  [2]  for  the  case 
of  longitudinal  impact,  without  torsion.  It  is  similar  to  the  theory  presented 
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by  Bleich  and  Nelson  [3]  for  the  case  of  a  uniformly  distributed  step  load  of 
pressure  and  shear  on  a  half  space,  except  that  Bleich  and  Nelson  assumed  the 
material  to  be  elastic-perfectly  plastic.  Rakhmatulin  [4]  and  Cristescu  [5]  have 
used  a  deformation  theory  of  plasticity  to  study  elastic-plastic  wave  propagation 
for  combined  stress  states  analogous  to  the  stress  states  considered  here.  They 
considered  the  case  of  edge  impact  of  two  plates  for  which  the  directions  of  the 
velocities  of  the  two  plates  are  in  the  plane  of  the  plates  and  are  oblique  to 
the  impacting  faces;  for  their  case  the  normal  strain  in  the  direction  parallel 
to  the  impacting  faces  and  in  the  plane  of  the  plates  is  zero  whereas  in  the  case 
considered  here  the  normal  stress  in  the  corresponding  direction  is  assumed  to 
be  zero.  Rakhmatulin  described  a  solution  in  which  first  a  uniaxial  stress  state 
propagates  into  the  body  as  a  simple  wave  followed  by  a  constant  state;  then  a 
combined  stress  state  propagates  as  a  wave  of  strong  discontinuity  followed  by  a 
constant  state.  Cristescu,  using  the  same  equations  as  Rakhmatulin  and  assuming 
a  material  for  which  the  slope  of  the  stress-strain  curve  in  shear  is  less  than 
the  secant  modulus,  found  two  types  of  combined  stress  waves.  He  showed  that  a 
given  level  of  compressive  strain  propagates  faster  in  the  faster  of  these  two 
combined  stress  waves  than  it  would  in  a  plastic  wave  for  a  uniaxial  stress  state; 

I  *  i  ' 

similarly 9  in  the  slower  combined  stress  wave  the  velocity  of  propagation  of  a 
given  level  of  shear  strain  was  shown  to  be  faster  than  in  a  plastic  wave  of  pure 
shear o  From  these  results  Cristescu  concluded  that  combined  dynamic  stress  is 
transmitted  in  a  body  only  by  combined  stress  waves.  In  contrast  with  the  results 
of  Rakhmatulin  and  Cristescu,  the  present  theory  predicts  that  a  wave  of . strong 
discontinuity  occurs  only  in  elastic  regions  and  that  a  leading  uniaxial  plastic 
wave  exists  for  ;a  step  load  of  normal  stress  and,  shear  *  A i bibliography  of  studies 
of  one-dimensional  plastic  waves  involving  combined  stress  states  is  given  in  a 

■  ,i  ! 

survey  paper  by  Cristescu  [6]. 
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1.  Derivation  of  the  Basic  Equations 

Consider  a  long  slender  thin-walled  tube  with  mean  radius  r  as  shown  in 

»  *  ‘ 

Fig.  1.  Let  U(x,t)  denote  the  average  displacement  in  the  longitudinal  direction 
at  time  t  of  the  cross-section  at  a  distance  x  from  the  impact  end f of  the  tube. 

Since  this  is  to  be  a  small  deformation,  theory,  the  coordinate  x.  can  refer  to 
either  the  initial  position  of  the  cross-section  or  its, position  at  time  t.  Let  ^ 

e(x,t)  denote  the  longitudinal  strain  U^,  (subscripts  x  and  t  denote  partial 
differentiation  with  respect  to  x  and  t  respectively)  and  u(x,t)  the  longitudinal 
particle  velocity1 U^.  Then  if  U  is  twice  continuously  differentiable 

et  *  ux  U) 

Analogously,  let  Q(x,t)  denote  the  average  rotation  about  the  x-axis  at  time  t 
of  the  cross-section  at  x0  Let  y(x>t)  denote  the  shearing  strain  rQ^,  and  v(x,t) 
the  tangential  velocity  rfit®  Then,  neglecting  changes  in  r,  we  have 

Yt  =  vx  (2) 

Conservation  of  linear  momentum  gives 

0X  =  put  (3) 

where  o(x,t)  is  the  force  in  the  longitudinal  direction  per  unit  cross-sectional 
area  and  p  is  the  density.  Conservation  of  angular  momentum  gives,  again  neglecting 
changes  in  r, 

tx  =  pyt  (4) 

where  x(x,t)  is  the  force  in  the  tangential  direction  per  unit  cross-sectional  area. 

The  strain  rates  and  y  are  assumed  to  be  the, sum  of  an  elastic  part 
and  a  plastic  part. 
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Thus , 


(5a) 


y 


t 


(5b) 


where  the  superscripts  e  and  p  denote  the  elastic  and  plastic  parts  respectively. 
In  formulating  the  stress-strain  relations  it  will  be  assumed  that  a  and  t  are 
the  only  non-zero  stresses;  in  particular,  the  hoop  stress  associated  with  the 
lateral  inertia  of  the  tube  will  be  neglected.  For  an  isotropic  elastic  solid 
the  elastic  parts  in  Eqs.  (5)  are  given  by 


1 

E  °t 


Y+-  = 


1 

M  Tt 


(6a) 


(6b) 


where  E  is  Young's  modulus  and  y  is  the  modulus  of  rigidity. 
The  plastic  strain  rates  are  given  by 


(7a) 


(7b) 


where  f(a,t)  is  the  yield  function  and  X  is  a  positive  scalar  function  to  be 
specified  subsequently.  The  yield  condition  is  assumed  to  have  the  form 

f(o,t)  =  (j)  +  t2  =  k2  (8) 

where  6  is  a  constant,  and  k  is  the  yield  stress  in  pure  shear.  For  6  =  /3  Eq.  (8) 
corresponds  to  the  von  Mises  (distortion  energy)  yield  condition  and  for  6=2 
Eq.  (8)  corresponds  to  the  Tresca  (maximum  shear  stress)  yield  condition.  If 
isotropic  work-hardening  is  assumed  then  there  is  a  one  to  one  correspondence 
between  k  and  the  plastic  work  W^*.  The  plastic  work  pate  is  given  by 
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+  xyp 


(9) 


From  Eqs.  (7),  (8),  and  (9)  and  the  relation  WP 
to  satisfy 


dw  •  . 

-rz —  k  the  function  X  is  found 
dk  t 


\  dWP 
2k2  dk 


(10) 


The  function  WP(k)  can  be  determined  from  the  stress-strain  curve  for  a 

'  1  ! 

simple  tension  test  or  a  pure  shear  test.  For  this,  let  o  =  f(e)  be  the  stress- 
strain  curve  in  simple  tension.  An  increment  in  stress  da  is  related  to  the 

g  P 

elastic  strain  increment  de  and  the  plastic  strain  increment  der  by 

da  =  f  (e)[dee  +  dep]  (11) 


Substituting  de6  =  da/E  and  dep  =  dWp/a  in  Eq.  (11)  and  regarding  the  slope  of  the 

t 

stress-strain  curve  (i.e.  f  (e))  as  a  function  of  the  stress  a,  say  g(a),  we 
obtain 


,  /  \  pda  dW  -i 

da  =  g  ( a )  [—  +  —  ] 


(12) 


For  simple  tension  the  yield  condition,  Eq.  (8),  reduces  to  a  =  0k.  Using  the 
latter  relation  in  Eq.  (12)  we  obtain  the  following  first  order  differential 
equation  for  the  function  WP(k). 


dWP 

dk 


02k  ( 


1 

g(0k) 


(13) 


Eq.  (13)  and  the  initial  condition  WP(Oy/0)  =  0,  where  is  the  initial  yield 
stress  in  tension,  determine  the  function  Wp(k)  corresponding  to  the  stress-strain 
curve  a  =  f(e). 

The  governing  equations  are  obtained  by  using  Eqs.  (5),  (6),  (7),  (8),  and 
(10)  to  express  and  y^  in  Eqs.  (1)  and  (2)  in  terms  of  a,  x,  a^  and  Tt.  Thus, 


6 


Eqs.  (1)  and  (2)  become 


|f  at  +  H(k)[(~)  at  +  ct  Tt]  =  ux 


—  x  +  H(k)[ax  a  +  02  t2  t  ],  =  v 
]}  t  t  t  x 


respectively ,  where 


If  Eq.  (8)  is  used  to  express  k  in  Eq.  (16)  in  terms  of  a  and  x  then  Eqs.  (3), 
(4),  (14),  and  (15)  constitute  a  system  of  four  first  order  partial  differential 
equations  for  the  four  unknown  functions  a,  t,  u  and  v* 


2.  Characteristic  Properties  of  the  Equations 


The  governing  system  of  first  order  partial  differential  equations  derived 
in  the  previous  section  can  be  written  in  the  following  matrix  form 


L[w]  =  Aw  +  Bw  =0 

L  X 


where  w  is  the  vector 


H  =  0 


and  A  and  B  are  the  following  symmetric  matrices , 


0  £+H(J) 


0  -  +  H62t2 

y 


o  -1 
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The  system  of  equations,  Eqs.  (17),  is  a  system  of  quasi-linear ,  symmetric  hyper¬ 
bolic,  partial  differential  equations  of  first  order.  (The  theory  of  such  equations 

i 

is  given,  for  example,  in  Chapter  5  of  Ref.  [7]).  The  characteristic  velocities  c 
for  Eqs.  (17)  are  the  roots  of  the  characteristic  equation  det.  C=0  where  C  =  (cA-B) 
is  the  characteristic  matrix.  The  characteristic  equation  is 

L(pc2 )  -  (M+N)(pc2)  +  1  =  0  (18) 

where 

L  =  ^  +  ^H(t,/e)2  + 1 1,02,2 ' 

M  =  —  +  H02t2 

y 

N  =  g-  +  H(0/0)2 
The  roots  of  Eq.  (18)  are 

pc 

The  wave  speed,  taking  the  plus  sign  inEq.  (19),  will  be  denoted  by  c^  and, 

taking  the  minus  sign,  by  c  ;  the  subscripts  f  and  s  denote  'fast'  and  'slow' 

s 

wave  speeds  respectively.  For  the  elastic  case  (i.e.  H=0)  c^  becomes  the  elastic 

1/2 

bar  velocity  cq  =  (E/p)  and  cg  becomes  the  elastic  shear  wave  velocity 
1/2 

c^  =  (y/p)  •  For  the  case  of  longitudinal  impact  (i.e.  x=0)  c^  corresponds 

to  the  velocity  of  longitudinal  plastic  waves  for  stresses  sufficiently  small 

that  c_  >  c9;  otherwise,  c  corresponds  to  the  velocity,  of  longitudinal  plastic 
r  s 

waves.  For  the  case  of  torsional  impact  (i.e.  0=0)  c  corresponds  to  the  velocity 

s 

of  torsional  plastic  waves. 

In  order  to  study  the  dependence  of  the  wave  speeds  on  the  stress  state 

■  ■  j 

it  is  helpful  to  characterize  the  stress  state  (o,t)  by  the  associated  value 
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of  the  yield  stress  k,  from  Eq.  (8),  and  an  angle  <j>  by  use  of  the  relations 


t  =  k  sin  <f> 

(20a) 

a/0  =  k  cos  <J> 

(20b) 

for  which  the  yield  condition,  Eq.  (8),  is  identically  satisfied.  Substituting 
Eqs.  (20)  in  Eq.  (19)  we  obtain  expressions  for  the  wave  speeds  in  terms  of  k 
and  These  expressions  can  be  simplified  by  substituting  for  H,  using  Eqs.  (16) 
and  (13),  and  writing  the  slope  of  the  uniaxial  stress-strain  curve  at  a  =  0k 
(i.e.  g(6k))  as  8E  where  6  =  6(k)  satisfies  0  <  8  <  1  and  decreases  monotonically 
with  increasing  k.  In  this  .way.  we  obtain  the  following  dimensionless  form  for  the 
wave  speeds  as  functions  of  8  and  <f>. 

1/2 

(_c_)2  _  M  t  N  t  [(M-N)2  +  P]  (21) 

C2  ;  :  2  L 

where 

2 

L  =  1  +  (1/ 8-1) (cos2  <j)  +  2 ( lVv)  Sin2  ^ 

M  =  2(l+v)  +  (1/8-1)  02  sin2  <j> 

N  =  1  +  (1/8-1)  cos2  <j> 

P  =  4(l/8-l)2  02  sin2  <j>  cos2  <f> 

and  v  is  Poisson's  ratio.  The  wave  speeds  c_  and  c  given  by  Eq.  (21)  are  plotted 

r  s 

vs.  <}>  in  Figo  2  for  fixed  values  of  0,  v,  and  $. 

For  values  of  v  satisfying  0  <  v  £  1/2  and  values  of  0  satisfying  3  £  02  <  4, 
the  wave  speeds  c  and  c  satisfy  c  £.  c9  £  c.  £  c  •  Furthermore 5  for  an  arbitrary 

t  S  S  a  I*  O 

fixed  value  of  <J>9  c  and  c^.  decrease  monotonically  with  decreasing  3.  Also  ?  for  a 

S  x 

fixed  value  of  8,  c  is  monotonically  decreasing  and  cf  is  monotonically  increasing 

s  r 

with  increasing  <j)  for  <j>  in  the  interval  (0,  tt/2). 
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The  characteristic  curves,  or  simpiy  the  characteristics,  for  Eqs.  (17) 
are  the  four  families  of  curves  which  satisfy 


dx  _ 

dt  "  *  cf 


(22a) 


dt  *  Cs 


(22b) 


where  cf  and  c  are  given  by  Eq.  (19).  Since  the  right  sides  of  Eqs.  (22)  depend 
r  s 

on  a  and  t  the  characteristics  can  be  constructed  only  after  the  solution  a(x,t) 
and  x(x,t)  is  known. 

Associated  with  the  characteristic  matrix  C,  for  a  characteristic  with  wave 
velocity  c,  is  a  null  vector  £  which  satisfies 


C£  =  0 


There  are  four  linearly  independent  null  vectors  corresponding  to  the  four  wave 
velocities  ±  c^  and  ±  cg.  The  null  vectors  for  Eqs.  (17)  which  correspond  to 

;  i 

positive  wave  velocities  are  given  by 


-pc^  Hat 


-i-  -  cJ±  +  H(a/0)2) 


-  dx 

for  dt  =  +  Cf 


(24a) 


pc|  (|  +  H(a/0)2)  -  1 


— - c  (i  +  H02t2) 


pc2  (-  +  H02t2)  -  1 
s  y 


dx 

for  TT  =  +  c 
dt  s 


(24b) 


-pc2  Hot 
s 
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in  which  there  is  an  arbitrary  scalar  factor.  The  null  vectors  and  i 

r  s 

corresponding  to  negative  wave  velocities  are  obtained  from  Eqs.  (24a)  and  (24b) 

respectively  by  replacing  c,.  by  -c,.  and  c  by  -c  . 

i  t  .  s  s 

The  functions  a,  x,  u,  and  v  are  continuous  across  a  characteristic,  but 
jumps  in  their  derivatives  may  occur  in  the  direction  of  the  normal  to  the  charac¬ 
teristic,  The  jump  [w^]  in  the  normal  derivative  of  the  vector  w  must  satisfy 
C[wn]  =  0  so  that,  in  view  of  Eq.  (23),  [wn]  is  proportional  to  the  corresponding 
null  vector  Since  all  components  of  the  null  vector  are  non-zero,  except  for 
special  cases  in  which  either  H=  0,  tr=  0,  or  t  =  0,  it  follows  that  in  general 
both  longitudinal  and  torsional  motion  are  associated  with  each  wave  speed. 

The  linear  combination  of  Eqs.  (17)  given  by 

£;L[w]  =0  (25) 

where  the  dot  denotes  the  Euclidean  inner  product,  is  an  interior  differential 
equation  on  the  characteristic  associated  with  the  null  vector  l  (Ref.  [7]  pp.  424- 
427).  Substituting  Eqs.  (24)  in  Eq.  (25)  we  obtain  the  following  incremental 
relations  along  the  characteristics. 

[1/pc2  -  (1/E  +  H(a/6)2)]  [pcf  dv  -  dx]  + 

Y 

HaxCpc,.  du  -  da]  =  0;  -r—  =  +  c  (26a) 

f  at  t 

[l/pc2  -  (1/p  +  H02x2)]  [pc  du  -  da]  + 
s  s 

A  v 

Hax[pc  dv  -  dx]  -  0;  -r—  =  +  c  (26b) 

s  ’  dt  s 

dx  dx 

The  incremental  relations  along  the  characteristics  =  -  c^  and  ~r r  =  -  c  are 

&  dt  f  ,  dt  s 

obtained  by  replacing  c^.  by  -c^  and  c  by  -c  in  Eqs.  (26a)  and  (26b)  respectively. 

X*  X  s  s 

These  incremental  relations  along  characteristics  are  used  in  the  next  section  to 
obtain  simple  wave  solutions.  They  are  also  useful  in  deriving  difference  equations 
for  the  numerical  solution  of  problems  involving  arbitrary  initial  and  boundary 
conditions  [8]. 
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3.  Simple  Wave  Solutions 


There  are  two  types  of  simple  wave  solutions  of  Eqs.  (17);  namely,  slow 

o  •  e  •  •  dx 

simple  waves  in  which  the  solution  is  constant  along  the  characteristics  =  cg5 

and  fast  simple  waves  in  which  the  solution  is  constant  along  the  characteristics 

dx  •  •  •  • 

—  =  c^.  Characteristics  for  these  two  types  of  simple  waves  are  shown  in  Fig.  3. 

The  solution  for  slow  simple  waves  can  be  obtained  by  examining  the  com¬ 
putation  of  the  solution  at  point  P  in  Fig.  3(a)  for  given  constant  values  of 

a,  t,  u9  and  v  on  the  characteristic  c°.  The  solution  at  point  P  is  computed 

s 

dx 

by  integration  of  Eqs.  (26)  along  the  characteristics  =  +  c^,  q.,c_p,  and  -  cg. 
Since  cr,  t,  u,  and  v  are  constant  along  c°  these  three  incremental  relations  can 
be  regarded  as  three  ordinary  differential  equations.  Adding  the  two  equations 
for  the  characteristics  +  c^  and  -  c^  gives  the  following  first  order  ordinary 
differential  equation  for  a  as  a  function  of  t 


't'(cr,T) 


where 


^(o,!) 


1/E  +  H(0/0)2  -  1/pc2 


(28a) 


or,  in  view  of  Eq.  ( 18 ),  equivalently 


4'(a,'r)  = 


1/p  +  H02t2  -  1/pc2 


(28b) 


It  is  convenient  to  also  have  expressions  for  ¥  involving  the  wave  speed  c  instead 

s 


of  c_p.  Thus 


T(a,t  ) 


1/y  +  H02t2  -  1/pc2 
s 


(28c) 
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or 

Hot 

’Kct.t)  = - 

1/E  +  H(a/6)2  -  1/pC2 

s 


( 28d ) 


From  Eq.  (28b),  V  vanishes  on  the  x-axis,  on  the  initial  yield  surface,  and  on  the 
a-axis  in  the  interval  0y  <  |a|  <  a  where  a  is  the  stress  for  which  the  slope  of 
the  uniaxial  stress-strain  curve  is  equal  to  y,  From  Eq,  (28a),  V  is  infinite  at 
points  on  the  a-axis  for  which  |a|  >  a.  Except  at  these  points,  the  function 
T(c,x)  is  continuous  and  continuously  differentiable Thus,  a  unique  solution 
of  Eq,  (27)  passes  through  any  point  (a,x),  except  the  points  where  V  is  infinite. 
This  solution  will  be  denoted  by  a  =  a(x;  (a,x)).  The  solution  through  a  point 
(a,x)  with  o  ?  0  cannot  intersect  the  x-axis  because  a  =  0  is  a  unique  solution  of 
Eq.  (27).  Likewise,  the  solution  through  a  point  (a,x)  with  x  t  0  cannot  inter¬ 
sect  the  a-axis  at  values  of  a  satisfying  |a|  >  a  because,  at  such  points,  x  =  0 

is  a  unique  solution  of  the  differential  equation  dx/da  =  l/¥(a,x).  Therefore, 
the  solution  of  Eq.  (27)  through  an  arbitrary  point  (a,x)  must  either  intersect 
the  initial  yield  surface  or  intersect  the  a-axis  in  the  interval  0y  £  |a|  <  a. 

Solutions  of  Eq.  (27),  obtained  by  numerical  integration,  are  shown  as  solid 
lines  in  Fig.  40  For  these  curves  a  stress-strain  curve  representing  that  of 
annealed  commercially  pure  aluminum  was  used.  The  von  Mises  yield  condition 
(i.e.  e  =  /3)  was  used. 

In  Fig,  4  there  are  three  types  of  curves:  (a)  curves  that  pass  through 

the  point  a  =  a,  x  =  0;  (b)  curves  that  intersect  the  a-axis  in  the  interval 

ay  <  | a |  <  a;  (c)  curves  that  intersect  the  initial  yield  surface.  Each  of  the 
curves,  or  stress  paths,  represents  a  locus  of  stress  states  that  occur  at  a 
section  during  the  passage  of  a  slow  simple  wave.  These  stress  paths  are  tra¬ 
versed  in  the  directions  indicated  by  the  arrows  in  Fig.  4.  The  slopes  of  the 
curves  in  Fig*  4  satisfy  the  inequality 


0  <  3—  <  — — — 

dT  e2i 


(29) 


and  approach  the  upper  limit  as  the  radius  of  the  subsequent  yield  surface, 
[(a/0)2  +  x2f  ,  becomes  infinite. 

The  solution  for  stresses  as  functions  of  x  and  t  for  a  slow  simple  wave., 
centered  at  the  origin,  is  determined  implicitly  by  the  relations 


a  =  a(x;  (0 ,t ) ) 


(30a) 


Cs(a’T)  =  t 


(30b) 


where  a  =  a(x;  (a,x))  is  the  stress  path  for  the  slow  simple  wave,  and  the  function 

c  (a,x)  is  given  by  Eq.  (19).  In  order  for  this  slow  simple  wave  solution  to  be 
s 

valid,  however,  the  wave  velocity  c  must  decrease  monotonically  as  the  stresses 

s 

a  and  x  increase  along  the  curve  a  =  a(x;  (a,x)).  To  show  that  this  is  the  case 
it  is  convenient  to  regard  cg  as  a  function  of  g  and  <$>  as  in  Eq.  (21).  Then 


_ _ s  _  _ s  c[g  dk  _ £  d$ 

dx  "  9g  dk  dx  9(f)  dx 


s  s 

Recall  that  -g-g—  is  positive  and  is  negative  (see  Fig.  2).  From  the  condition 
that  the  stress-strain  curve  is  concave  downward  d3/dk  is  negative.  The  result 
that  do/dt  satisfies  inequality  (29)  insures  that  dk/dx  and  d<f>/dx  are  positive. 
Consequently,  dcg/dx  is  negative  and  the  simple  wave  solution  is  valid. 

The  slow  simple  wave  solution  is  completed  by  determining  the  velocities 
u  and  v  associated  with  each  point  on  the  stress  path  a  =  a(x;  (a,x)).  For  this, 
use  is  made  of  the  two  ordinary  differential  equations  which,  in  addition  to 


Eq.  (27),  result  from  integrating  Eqs.  (26)  along  the  characteristics 
-  c^,  and  -  cg  in  Fig.  (3a).  These  two  equations  can  be  written  as 


=  +  c-. 


(32a) 


du 

da 


(32b) 
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Since  the  right  side  of  Eqs.  (32)  is  a  known  function  of  0  and  t  the  change  in 
u  and  v  as  the  stress  state  varies  along  the  path  a  =  a(x;  (a,x))  can  be  obtained 
by  quadrature,, 

The  change  in  strain  with  the  passage  of  a  slow  simple  wave  can  be  obtained 
by  integration  of  the  strain  rates  *  Eqs,  (14)  and  (15),  with  respect  to  time, 
along  the  path  a  =  a(x;  (a,x)).  Making  use  of  Eqs.  (27)  and  (28),  we  can  write 
the  equations  governing  the  change  in  strain  in  the  following  form 


de  1 


dy 

dt 


(33a) 

(33b) 


from  which  the  strains  can  be  obtained  by  quadrature.  An  interesting  consequence 
of  Eqs.  (33)  is  that  in  slow  simple  wave  regions 


de  _  d£ 
dy  dx 


(34) 


That  is,  the  direction  of  the  total  strain-rate  vector  coincides  with  that  of  the 
stress-rate  vector. 

The  solution  for  fast  simple  waves  is  obtained  in  the  same  way  as  was  the 
solution  for  slow  simple  waves.  Integration  alopg  the  three  characteristics  in 
Fig.  3b  yields  the  following  differential  equation  for  the  stress  path  for  fast 
simple  waves 


da  _  1 

dx  ~  yCa,!) 


(35) 


where  ¥(a,x)  is  the  same  as  in  Eq.  (28).  Comparison  of  Eqs.  (35)  and  (27)  reveals 
that  the  stress  paths  for  fast  simple  waves  are  orthogonal  to  the  stress  paths 
for  slow  simple  waves.  The  former  are  shown  as  dashed  lines  in  Fig.  4.  As  these 
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stress  paths  are  traversed  from  left  to  right  the  wave  speed  decreases  mono- 
tonically  from  cq  to  c2.  The  equations  determining  the  changes  in  velocities  for 
fast  simple  waves,  analogous  to  Eqs.  (32),  are 


du  _  1 

da  "  pcf 

dv  _  _  1 

dx  -  pc^ 


The  equations  for  computing  the  changes  in  strains  are 


(36a) 


(36b) 


de 

da 


(37a) 


dy 

dx 


(37b) 


From  Eqs.  ('37),  the  directions  of  the  total  strain-rate  vector  and  the  stress-rate 
vector  also  coincide  for  fast  simple  waves. 


4.  Solution  for  Step-Loading  of  a  Semi-Infinite  Tube 


Solutions  for  a  variety  of  step-loading  problems  for  a  semi-infinite  tube 
can  be  obtained  by  appropriate  combinations  of  slow  and  fast  simple  wave  solu¬ 
tions,  constant  state  solutions,  and  waves  of  strong  discontinuity  (i.e.  waves 
across  which  jumps  in  stress  and  velocity  occur).  Jumps  in  a  and  u  can  occur  only 
across  a  wave  front  with  speed  cq;  jumps  in  x  and  v  can  occur  only  across  a  wave 
front  with  speed  c2„  From  conservation  of  momentum  these  jumps  must  satisfy 


[a]  =  -  pc  [u]  for  =  c 
o  dt  o 


[t]  =  -  pc2Cv]  for  =  c2 


(38a) 


(38b) 


16 


where  square  brackets  denote  the  jump  in  the  enclosed  quantity. 

As  an  example,  we  consider  the  case  of  a  tube  initially  at  rest  and  unstressed 

which,  at  the  end  x  =  0,  is  simultaneously  subjected  to  a  constant  normal  stress 

a  and  a  constant  shear  stress  t  „  That  is,  we  seek  a  solution  of  Eqs.  (17)  in 
o  o 

the  region  0  <  x  <  0<t<°°  which  satisfies  the  initial  conditions 

u(x,0)  =  a(x,0)  =  v(x,0)  =  t(x,0)  =  0  (39) 

and  the  boundary  conditions 


a(0,t)  =  a  0  <  t  <  <=° 

o 

(40a) 

T(0,t)  =  T  0  <  t  <  ” 

(40b) 

o 

in  which  and  are  constants  satisfying 

(o  /e)2  +  t2  .>  (a  /e)2 

o  o  Y 

(41) 

where  a y  is  the  initial  yield  stress  in  simple  tension. 

For  this  problem  there  are  three  types  of  wave  solutions  as  shown  in  Figs. 

5(a),  5(b),  and  5(c);  which  type  results  in  a  given  case  depends  on  the  values 
of  and  t  .  Figs.  5(a),  5(b),  and  5(c)  correspond  to  the  cases  when  a ^  and 
are  such  that  the  point  (aQ,T  )  in  Fig.  4  is  intersected  by  solutions  of  Eq.  (27) 

of  types  (a),  (b),  and  (c),  respectively.  Double  parallel  lines  in  Fig.  5  repre¬ 

sent  wave  fronts  across  which  jumps  in  stress  and  velocity  occur.  For  case  (a) 
the  stress  wave  behavior  is  as  follows:  a  jumps  to  a y  across  the  wave  front  pro- 

a 

pagating  with  a  speed  cq;  a  increases  across  the  fast  simple  wave  region  to  a, 
according  to  the  von  Kcirman  solution,  while  T  remains  zero;  a  and  t  increase 
across  the  slow  simple  wave  region,  following  stress  path  (a)  in  Fig,  4,  until 
the  stress  state  o  =  a  ,  t  =  t  is  reached;  a  and  t  are  constant  in  the  constant 
state  region.  Case  (b)  differs  from  case  (a)  only  in  that  the  fast  simple  wave 
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solution  terminates  when  the  stress  o  reaches  the  stress,  say  at  which  the  stress 
path  (b)  in  Fig*  4  intersects  the  0-axis.  A  constant  state  region  with  a  = 
and  t  =  0  separates  the  fast  and  slow  simple  waves  in  Fig*  5(b)*  In  case  (c)  a 
jumps  to  (see  Fig.  4)  across  the  wave  front  x  =  CQt ;  t  jumps  from  0  to  across 

the  wave  front  x  =  c^t;  a  and  t  increase  across  the  slow  simple  wave  region, 

■i 

following  stress  path  (c),  to  the  stress  state  0  =  oo>  t  =  t  .  The  particle 

l 

velocities  u  and  v  for  these  three  causes  can  be  obtained  from  the  solution  for 

stresses  by  use  of  Eqs.  (32),  (36),  and  (38)* 

As  a  second  example  we. consider  a  case  which  differs  from  the  first  example 

only  in  that  the  tube  is  initially  stressed  statically  by  a  shear  stress  t  which 

exceeds  the  initial  proportional  limit  of  the  material*  In  this  example  the  final 

shear  stress  will  again  be  denoted  by  tq  so  the  jump  in  shear  stress  at  t  =  0  is 

t  -To  The  solution  for  this  case  can  be  read  from  Figs.  6(a)  and  6(b);  a  in- 
o  s 

creases  and  t  decreases  across  the  fast  simple  wave  region,  following  the  stress 
path  ab,  to  the  stress  state  a  =  <S^9  x  =  there  is  a  constant  state  region  with 
a  =  t  =  for  cs(ai,T1)  <  ^  <  c^(a1>x1);  a  and  t  follow  the  stress  path  be 
as  they  increase  across  the  slow  simple  wave,  region  to  the  final  stress  state 

(oo>To)< 

This  second  example  exhibits  the  unexpected  result  that  the  shear  stress 
at  any  cross-section  would  first  decrease  due  to  the  passage  of  a  fast  simple  wave. 

It  should  be  possible  to  carry  out  experiments  that  would  either  verify  or  refute 
such  a  prediction* 

A  case  for  which  the  solution  exhibits  an  even  more  unexpected  behavior  is 
the  case  where  the  tube  is  initially  stressed  by  a  tensile  stress  Gg  and  then 
subjected  to  step-loading  such  that  the  final  stress  state  (0^,1^)  lies  outside 
the  yield  surface  through  the  point  (0^,0).  That  is,  o and  satisfy 

(a  /e)2  +  t2  >  Co  /e)2 
o  o  s 


(42) 
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The  solution  for  this  case  is  shown  in  Fig.  7.  Unloading  occurs  at  the  wave 

front  x  =  CQt  where  the  normal  stress  falls  to  ;  x  jumps  to  x ^  across  the  wave 

front  x  =  c2t;  a  and  x  follow  the  stress  path  cd  as  they  increase  across  the  slow 

simple  wave  region  to  the  final  stress  state  (aQ, xq).  In  this  case  the  magnitude 

of  the  jumps  at  the  two  wave  fronts  clearly  depends  strongly  on  the  assumption 

that  the  subsequent  yield  surface  for  tension  followed  by  torsion  can  be  located 

by  the  theory  of  isotropic  work-hardening.  However,  in  static  tests  the  subsequent 

yield  surface  ordinarily  agrees  better  with  the  assumption  that  the  yield  surface 

translates  in  the  direction  of  the  initial  tensile  stress  than  it  does  with  the 

assumption  of  isotropic  work-hardening  [9],  Therefore,  the  theoretical  predictions 

in  this  case  probably  do  not  agree  well  with  what  would  be  observed  experimentally. 

So  far  solutions  have  been  given  for  step-loading  of  tubes  subjected  to 

either  an  initial  shear  stress  x  =  x  or  an  initial  normal  stress  a  =  a  .  The 

s  ,  s 

solution  for  the  case  of  an  arbitrary  initial  stress  state  cr  =  a  ,  x  =  x  can 

S'  s 

clearly  be  obtained,  without  difficulty,  using  the  techniques  already  established 
here , 

The  solutions  given  here  are  also  valid  for  the  case  of  combined  longitu¬ 
dinal  and  torsional  impact  in  which  the  velocities,  instead  of  the  stresses,  are 
prescribed  at  x  =  0.  For  given  values  of  prescribed  velocities  (u^v^)  the  cor¬ 
responding  stresses  at  the  impaet^face  (0^,1^)  can  he  obtained  by  interpolation 

from  a  network  of  contours  of  constant  values  of  u  and  v  as  shown  in  Figo  80  In 

o  o 

Figo  8a  contour  u^  ~  constant  represents  the  locus  of  applied  stresses  (0^,1^) 
for  which  the  longitudinal  velocity  at  the  impact  end  is  u^.  The  contours  v^  =  con¬ 
stant  have  an  analogous  interpretation o  These  contours  are  obtained  by  first 
determining  several  stress  paths  for  slow  simple  waves  as  shown  in  Fig6  4  and  then 
computing  the  particle  velocities  along  these  paths  by  numerical  integration  of 
Eqs<>  (32)0  The  contours  are  then  plotted  by  connecting  points  on  the  various  stress 
paths  for  which  the  velocities  are  the  same0 
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FIG.  5  CHARACTERISTICS  FOR  STEP -LOADING 


FIG.  6  SOLUTION  FOR  STEP-LOADING  OF  A  PRE-TWISTED  TUBE 
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FIG.  7  SOLUTION  FOR  STEP-LOADING  OF  A  PRE-TENSIONED  TUBE 


ASSUMED  MATERIAL  PROPERTIES  THE  SAME  AS  FOR  FIG. 4 
,  =  CONSTANT  (  IN./SEC.)  V0  =  CONSTAN  T(  IN./SEC.) 


